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Quantitative density measurements from electron scattering show that shear bands in deformed 
AlgsYyFes metallic glass exhibit alternating high and low density regions, ranging from -9 % to +6 % 
relative to the undeformed matrix. Small deflections of the shear band from the main propagation 
direction coincide with switches in density from higher to lower than the matrix and vice versa, 
indicating that faster and slower motion (stick slip) occurs during the propagation. Nanobeam 
diffraction analyses provide clear evidence that the density changes are accompanied by structural 


changes suggesting that shear alters the packing 
clusters. This bears a striking resemblance to the 
upon deformation of granular media. 

PACS numbers: 61.05.J-, 61.43.Dq, 83.50 

Metallic glasses (MGs) offer unique properties such as 
high strength, extended elasticity, high wear and corro¬ 
sion resistance, or excellent soft magnetic behavior [1]. 
However, the limited ductility and especially the immedi¬ 
ate catastrophic failure in tension once the elastic limit is 
reached are major obstacles to applications as structural 
materials [1]. This behavior has led to substantial effort 
towards understanding and improving the accommoda¬ 
tion of plastic deformation in MGs, including the design 
of composites consisting of ductile crystalline phases in 
a bulk MG matrix [2-6]. Monolithic glasses with higher 
Poisson’s ratios near the ideal value of 0.5 [7-9] tend to 
have improved compressive and bending plasticity, and 
are reported to show a fine dispersion of shear bands [10]. 
Well below the glass transition temperature, plastic flow 
in metallic glasses is restricted to narrow regions called 
shear bands [11-13], widely believed to be dilated zones 
of increased free volume caused by shear localization en¬ 
abling shear softening [14-16]. Thus, the density of or 
free volume inside shear bands is often treated as key 
to understanding plasticity in MGs. Quite a large range 
of density or free volume changes have been observed, 
based on various experiments [15-23]. However, most 
experimental methods, including positron annihilation 
spectroscopy and calorimetry, measure the free volume 
integrated over the entire sample to be typically from 
+1 % to +3 %. High resolution transmission electron 
microscopy has shown the presence of nanovoids inside 
shear bands, believed to form from the coalescence of 
atomic-sized free volumes after deformation [24], but so 
far has not been used to measure material density either 
on average or locally within the shear band. 

Previously the local density within shear bands in 
AlggYyFes has been probed at selected positions and for 
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different shear bands using quantitative HAADF-STEM 
on deformed and undeformed metallic glass regions (see 
Supplemental Material for details). These experiments 
showed both positive and negative density changes for 
shear bands with respect to the undeformed glass [25]. 
Minor compositional changes relative to the matrix were 
observed in the shear band segments. While composi¬ 
tional changes contributed mainly to the positive density 
variations in the bright shear band segments, the change 
in free volume was the dominant effect for the dark shear 
band segments. 

Here we show, by combining several quantitative elec¬ 
tron scattering signals, that shear bands in AIggYyFes 
metallic glass have segments of both decreased and in¬ 
creased density that alternate along the propagation di¬ 
rection. The changes in density are correlated to small 
deflections in the propagation direction and changes in 
medium-range structural order and chemical composi¬ 
tion. Similar behavior has been reported for granular 
media, but not for amorphous solids. One important im¬ 
plication is that individual shear bands apparently prop¬ 
agate via local, segment-wise stick slip that might origi¬ 
nate from their complex topology. These results are im¬ 
portant to the physics of deformation in metallic glasses, 
and suggest connections to the physics of granular mate¬ 
rials and jammed systems. 

Ingots of the target composition AIggYyFes were pre¬ 
pared by arc melting pure A1 (99.999 %), Fe (99.99 %) 
and Y (99.9 %). Fully amorphous ribbons with a thick¬ 
ness of about 40 fi m were produced by melt spinning 
using a tangential wheel (Cu) speed of 47 ms -1 . The 
melt-spun ribbons were cold-rolled to a thickness reduc¬ 
tion of 23 % and subsequently prepared for transmis¬ 
sion electron microscopy (TEM) by twin-jet electropol- 
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ishing using HN 03 iCH 3 0 H in a ratio 1:2 at 253 K. The 
TEM study was performed with a FEI Titan 80-300 
(probe) aberration-corrected transmission electron mi¬ 
croscope operated at 300 kV in the scanning transmission 
electron microscopy (STEM) mode. The following condi¬ 
tions were used during the experiments: a probe current 
of 15 pA, a collection semiangle of 51.3 mrad for the high- 
angle annular dark-field (HAADF) detector (130 mm 
camera length), a convergence semiangle a of 0.6 mrad, 
a nominal spot size of 2 nm, and a step size of 2 nm. For 
electron energy loss spectroscopy (EELS) we used a col¬ 
lection semiangle /3 of 5.73 mrad, an entrance aperture of 
2.5 mm, an energy dispersion of 0.1 eV/channel and an 
acquisition time of 0.025 s. The fluctuation electron mi¬ 
croscopy (FEM) data were obtained from energy-filtered 
nanobeam diffraction using a camera length of 510 mm. 

HAADF-STEM images from deformed AlgsYyFes, 
such as shown in Fig. 1, frequently showed shear bands 
with contrast changes from bright to dark and vice versa 
along their propagation direction. Each contrast reversal 
is accompanied by a slight deflection [25], always within 
±5° of the main propagation direction of the shear band. 
The propagation direction is determined from the bifur¬ 
cation. The angular deflection range matches the angle 
between the shear bands and the direction of the applied 
shear stress reported in literature [13]. 

Figure 2 shows the foil thickness profiles from part of 
the sample area shown in Fig. 1 determined from low-loss 
EELS. There is a linear increase in foil thickness along 
the propagation direction of the shear band which is re¬ 
lated to the wedge shape of the TEM specimen. There is 
also a slight thickness variation between the two sides of 
the shear band, probably due to the shear offset. Thick¬ 
ness profiles across the shear band show that there is no 
preferential thinning of the shear band by TEM sample 
preparation. 

The density variation Ap of each bright or dark seg¬ 
ment was determined from the HAADF intensity and 
EELS thickness, as described in detail in the Supplemen¬ 
tal Material. The left side of Fig. 1 shows the local den¬ 
sity inside the shear band along its propagation direction. 
Repetitive density changes are observed. The density dif¬ 
ferences between the shear band and undeformed matrix 
range from (—9.0 ±1.0) % to (+6.3 ±1.0) %. The matrix 
region outside the shear bands does not show any impact 
of the deformation or of preexisting density fluctuations. 
The integrated density change along this specific shear 
band is A p = (—2.9 ± 4.3) %. No characteristic segment 
length was found over examination of many shear bands. 
However, the dark low density segments dominate the 
shear bands in general, which results in a negative value 
for the integrated density change. The mean density is 
in good agreement with previous reports of shear band 
density from macroscopic experiments [15, 22, 23]. The 
width of the shear band was measured on the clearest 
segments (//, ///, and IV) by fitting Gaussian peaks to 
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FIG. 1. Density variations in a shear band. Right: A 
HAADF-STEM image showing a shear band in cold-rolled 
AlssYrFes propagating from bottom to top with several con¬ 
trast changes. The individual segments are numbered (I — V). 
The deflection angles are indicated for the segments. Left: 
Quantified density variations of the shear band along the 
propagating direction with respect to the undeformed ma¬ 
trix. Bright and dark parts of the shear band are indicated 
by the white and gray background, respectively. 


the data in Fig. 1. The FWHM of wu = (3.72±0.21) nm, 
wm = (4.26 ± 0.16) nm and wiy = (4.24 ± 0.15) nm. 
Segment V is wider by inspection, so the consistency of 
the width in the other segments shows that they image 
the shear band near edge-on. When the shear band is 
inclined, the projected density difference along the elec¬ 
tron beam direction is suppressed due to the overlap of 
shear band and matrix. 

At this point it might be appropriate to question 
whether the shear band structures analyzed resemble the 
structures formed during the shear band propagation. 
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FIG. 2. Foil thickness calculations. Profiles of the av¬ 
eraged foil thicknesses for the boxed regions of the HAADF- 
STEM image displaying different shear band parts (see seg¬ 
ments //, ///, and IV in Fig. 1). The gray shaded areas 
indicate the position of the shear band. 


Since no crystal formation was found for the denser re¬ 
gions of the shear band, where crystallization may be 
expected to be facilitated [26], it seems safe to assume 
that the shear bands analyzed in this work do closely 
resemble the structures present during deformation. 

FEM provides information about the nanoscale or¬ 
der of the metallic glass from systematic coherent elec¬ 
tron nanodiffraction (see Supplemental Material) [27, 28]. 
Nanoscale order leads to variability in nanodiffraction 
into different directions in reciprocal space, which is cap¬ 
tured quantitatively in the annular mean of the variance 
image ^vimage(^) [27]. Figure 3 shows FEM data from 
a bright high-density segment (segment /), a dark low- 
density shear band segment (segment //), and the un¬ 
deformed matrix about 20 nm away from the dark shear 
band segment. The results for the lower density, dark seg¬ 
ment are similar to our previous report [25], but we have 
obtained better statistical quality data for the higher 
density, bright segment here. After accounting for dif¬ 
ferences in sample thickness, the variance signal for the 
matrix is consistent with previous FEM studies of the 
same glass composition [28], which data showed that un¬ 
deformed AlssYyFes contains nanoscale ordered regions 
with an internal structure similar to fee A1 embedded in a 



FIG. 3. Fluctuation Electron Microscopy. Annular mean 
of variance image Hyimage of different NBDP ensembles (un¬ 
deformed matrix, bright/dark parts of the shear band corre¬ 
sponding to the segments I and II). 


more disordered matrix. The primary peak in Hvimage(&) 
sits near the A1 { 111 } reflection, there is a shoulder at 
higher k consistent with A1 {200}, and there are no peaks 
at higher k. flyimage(^) for the bright, higher density 
shear band segment is similar to the matrix, but the 
{200} shoulder is suppressed. This indicates a similar 
structure of Al-like regions embedded in a matrix, but 
with reduced size or internal structural order for the Al- 
like regions [25, 28]. ^vimage(^) for the dark, low den¬ 
sity shear band segment indicates strikingly high struc¬ 
tural order. In addition to clear {111} and {200} peaks, 
there are peaks at the {311} and {220} A1 reflections as 
well [25]. However, all peaks are shifted towards lower 
k values. The reason for this has been shown to be hy¬ 
drostatic tensile strain originating from the surrounding 
amorphous material in which the Al-like regions or clus¬ 
ters are embedded [28]. The increase in Hvimage(^) shows 
that the Al-like regions are larger than in the matrix, 
and the persistence of peaks to high k shows that they 
have a higher degree of crystalline structural order [28]. 
Strong Al-like order is consistent with previous observa¬ 
tions of A1 crystallite nucleation in shear bands [26]. Here 
we emphasize that we have a mixture composed of Al- 
rich cluster/crystallites embedded in amorphous mate¬ 
rial present in the dark shear band segments rather than 
fully crystallized regions since the latter would produce 
enormous peaks in the fluctuation signal which we do 
not observe. Moreover, a former study showed that fully 
crystallized AlssYyFes metallic glass is composed of three 
crystalline phases, that is, fee Al, A1 3 Y (DO 19 structure) 
and a ternary intermetallic A^FesY phase [29]. Since the 
diffraction peaks of the FEM study closely match fee Al 
with the other two crystalline phases missing, a complete 
crystallization of the dark shear band segment can be ex¬ 
cluded here. We estimate the size of the ordered Al-like 
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regions to be <2 nm from the similarity of adjacent nan¬ 
odiffraction patterns [30]. Since we used a 2 nm diameter 
probe for the experiments, 2 nm is an upper bound and 
the real ordered regions may be even smaller. Moreover, 
for 3d transition metals such sizes are commonly associ¬ 
ated with clusters rather than with the crystalline state. 

Larger density decreases of 9 % or more as reported 
here must be accommodated by some local chemical 
changes, not just an increase in mean atomic spac¬ 
ing. A qualitative chemical analysis based on the Al-L 
edge showed an increase in the A1 signal for the dark 
shear band segments and a decrease for the bright ones 
[Fig. 4(d), Supplemental Material] which is consistent 
with our previous report [25]. It is worth noting that the 
chemical changes observed are minor. They are typically 
balanced between A1 and Fe. It was further shown that 
the A1 enrichment and the corresponding Fe depletion 
were not sufficient to account for the total density change 
in the dark shear band segments.The dominating effect 
causing the density change was the variation in free vol¬ 
ume. Based on the fact that Fe is not soluble in fee Al, we 
can estimate the maximum density variation arising from 
compositional changes. If we replace all of the present Fe 
(5 at%) by Al, it would account for a maximum density 
decrease of roughly 4.5 % which is not enough to explain 
a density decrease of 9 %. Thus, the density drop for the 
dark shear band segments must include changes in free 
volume in addition to the changes in composition [25]. 
Further information can be obtained from the zero los- 
sand plasmon signal extracted from the EELS data. The 
increased zero losssignal [Fig. 4(b)] for the dark shear 
band regions indicates less scattering events due to less 
dense material at a constant foil thickness. On the other 
hand, the plasmon peak [Fig. 4(c)] reveals an increase 
for such dark shear band segments, which also indicates 
that compositional changes (Al enrichment) must have 
occurred, since an Al increase enhances the free-electron 
density and thus the plasmon signal. 

The observation of an increased EELS ZLP for the dark 
shear band regions [see Supplemental Material Fig. 4(a)] 
and enhanced diffusion along shear bands by more than 
6 to 8 orders of magnitude compared to diffusion in the 
undeformed matrix [31] are perfectly consistent with en¬ 
hanced free volume in the shear bands. 

Structural changes inside shear bands have been ana¬ 
lyzed indirectly by low-temperature heat capacity mea¬ 
surements on severely deformed metallic glasses [32, 33]. 
The so-called Boson peak revealed that the atomic struc¬ 
ture was only modified inside the shear bands. Addi¬ 
tional studies of the aging dependence of the Boson peak 
indicated that two different regions exist in the deformed 
specimens. One of these regions contributes to an ac¬ 
celerated aging of states since after identical relaxation 
treatment between room temperature and the glass tran¬ 
sition, a state with lower enthalpy was attained by the 
deformed glass as compared to material without shear 


bands. Accelerated aging would be expected for the low- 
density regions of the shear band. 

The HAADF image intensity shows no evidence for the 
formation of voids, but at this sample thickness it is in¬ 
sensitive to the small (1 nm diameter) voids previously 
reported [24]. Reasons for void formation are the pres¬ 
ence of multiple deflections of the shear bands, stress con¬ 
centrators and low-viscosity regions. Nanovoids would 
also explain the early and catastrophic failure as well as 
the observation of vein patterns [8, 13] at the fracture 
surfaces of metallic glasses. 

The repeated switching between higher and lower den¬ 
sity observed along the shear band is directly analogous 
to the deformation behavior of granular media. Fazekas 
et al [34] have shown by simulations that deformation 
of idealized granular media composed of spheres led to 
the formation of shear bands with high and low density 
regions. They further showed that the critical density of 
the shear bands depends on friction which, in turn, de¬ 
fines a specific packing state. Thus the theoretical find¬ 
ings for granular media explain our observation of higher 
and lower density shear band regions as the result of fric¬ 
tion causing faster and slower motion of the shear bands. 
However, in the absence of granules, what plays their role 
in the metallic glass? It could be atoms, but the persis¬ 
tence of a similar structure of Al-like regions in a more 
disordered matrix within both the dark and bright seg¬ 
ments leads us to suggest that “the granules” may instead 
be tightly bound short- or medium-range atomic clusters, 
which can fill the space more or less efficiently [35-37]. 
Shear alters the packing of those clusters, but does not 
strongly modify their internal structure. 

Since glasses do not have a defined slip plane it is ex¬ 
pected that shear bands have a distribution of topologi¬ 
cal minima and maxima (“hills and valleys”) along their 
propagation direction, arising from inhomogeneity inher¬ 
ent to the glass. This nonplanar topology would give 
rise to a distribution of activation barriers for continued 
slip along the shear band. It is thus feasible, as a hy¬ 
pothesis, that compressive regions evolve before a local 
topological maximum and dilated regions may form after 
a topological maximum has been surmounted. Segments 
approaching a local barrier caused by a topological max¬ 
imum would propagate slower while the motion speed 
would increase again after passing such a barrier result¬ 
ing in a stick slip motion. This stick slip motion causes, in 
our opinion, the observed deflections. This interpretation 
fits very well with the observation in granular media [34] 
insomuch as the propagation velocity defines the packing 
state. 

In summary, we have quantified the density inside a 
shear band of deformed AlgsYyFes metallic glass along 
its propagation direction. The local density varies enor¬ 
mously, ranging from -9 % to +6 %, compared to the un¬ 
deformed matrix. Slight deflections from the main prop¬ 
agation direction of the shear band are found to coin- 



5 


cide with the density switching from positive to negative 
values. The density alterations are attributed to com¬ 
positional and free volume changes reflecting a different 
dense packing of atomic clusters, probably caused by dif¬ 
ferent propagation speeds. These findings bear a striking 
resemblance to granular shear bands formed upon defor¬ 
mation of granular media [34]. 
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SUPPLEMENTAL INFORMATION 


Density measurements 


We derive information about local density changes 
from the high-angle annular dark-field scanning transmis¬ 
sion electron microscopy (HAADF-STEM) signal (elec¬ 
trons collected by the HAADF detector) [25]. The dark- 
field intensity contains information about the mass 
thickness pt as follows: 


— = 1 - exp 

I o 


^ NaojA ^ 


1 — exp 



and for small arguments: 


I pt 
Io %k 


(1) 


where Na is the Avogadro’s number, a is the total scat¬ 
tering cross-section, p is the density, t is the foil thickness 
and A is the atomic weight. Xk is the contrast thickness, 
which is defined as . A simultaneously acquired EEL 
signal allows calculation of the specimen foil thickness t 
from the low-loss spectral region [38] and hence the deter¬ 
mination of density changes. Using Eq. (1) the relative 
density change (normalized to the matrix) may be writ¬ 
ten: 


A p 


Psb ~ Pm 
Pm 


IsBtM%k B 

Im^sb^ 


( 2 ) 


where pM , Psb are the mass densities, Jm, Isb are the 
HAADF intensities, x^f, x BB are the contrast thicknesses 
and tM,tsB are the corresponding foil thicknesses for the 
matrix and the shear band, respectively. The assumption 
of a constant contrast thickness Xk for matrix and shear 
band, which cause the Xk terms to cancel in Eq. 2, has 
been discussed in detail elsewhere [25]. The calculated 
foil thicknesses (Fig. 2) for the present investigation have 
been taken as input for calculating the densities in Eq. 
( 2 ). 


the collection semiangle, 6e is the characteristic scatter¬ 
ing angle of inelastic scattering corresponding to an en¬ 
ergy loss E , and Eq is the microscope voltage in kilovolt 
units. Calculations using the refractive index-corrected 
Kramers-Kronig sum rule are performed automatically, 
for example, within the Digital Micrograph software rou¬ 
tines. It analyses the single scattering distribution S ( E ), 
which is obtained from the EEL spectrum by removing 
the ZLP and plural inelastic scattering using the Fourier- 
Log method. For data processing the energy-shift of each 
individual spectrum has been corrected using DM scripts. 
The thickness at each point has been used in the density 
calculation (see Eq. (2)). 


Zero loss and Plasmon Peak 

The plasmon peak maximum was found at 15.2±0.1 eV 
by fitting the center of the zero loss peak (ZLP) and the 
plasmon peak (Lorentzian fits) and subsequently measur¬ 
ing their energy difference. This value confirms results 
from a FIB-prepared TEM sample [25]. An energy shift 
of the plasmon peak between the shear band segments 
and the matrix was not observed within the accuracy 
limit of ±0.1 eV (Fig. 4a). The individual profiles of the 
Zero loss and plasmon peak are shown in Fig 4b and 4c. 
Note that the dark shear band regions show an increase 
in the zero lossintensity which can only be explained by 
less dense matter. The opposite observation is found for 
the bright shear band regions which show decreased Zero 
loss intensities due to denser material. The same state¬ 
ments can be made for the plasmon peaks (Fig. 4c) which 
is consistent with our previous report [25]. A chemical 
analysis is added here (Fig. 4d) extracted from the A1 
L-edge using an energy window ranging from 73.5 - 81.5 
eV. We find an A1 enrichment for the dark shear bands 
where the pronounced fcc-like structures were observed 
in the FEM analysis (Fig. 3), whereas the bright shear 
band regions display a slight A1 decrease. 


Foil thickness calculation 

The refractive index-corrected Kramers-Kronig sum 
rule (Eq. (3)) was used to calculate the foil thickness 
based on the low-loss part of each individual EEL spec¬ 
tra [38]: 

4 ajm, I S(E)6E (3) 

(I-^)IzlpJ E]n(l + §P) 

where S (E) is the single scattering distribution, Izlp is 
the integrated intensity of the zero loss peak (ZLP), n 
is the refractive index of the material (~ 500 for met¬ 
als), F is a relativistic factor, ao is the Bohr radius, (3 is 


Influence of the collection semiangle on the HAADF 

signal 

The collection semiangle affects the collection of elec¬ 
trons contributing to the HAADF signal. This influence 
has been investigated along the typical wedge-shape of 
a TEM specimen. The results are shown in Fig. 5. It 
was found that for collection semiangles larger than 50 
mrad, the HAADF signal is essentially constant mean¬ 
ing that the majority of Rutherford-scattered electrons 
are properly detected. In the present analysis a collec¬ 
tion semiangle of 51.3 mrad (130 mm camera length) was 
used because it offered the best signal to noise ratio. 
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a) 



FIG. 4. Zero loss and plasmon peak/map (a) EEL 

spectra for different regions (undeformed matrix, bright/dark 
shear band segments; see nomenclature in Fig. 1). The gray 
regions indicate the size of the energy window used for ex¬ 
tracting the maps, (b), (c) Corresponding Zero loss and plas¬ 
mon profiles for the different regions. Right: Corresponding 
regions (boxes) in the HAADF-STEM image, (d) Profiles of 
the Al-L edge (energy window ranging from 73.5 - 81.5 eV) 
extracted from the individual EEL spectra. Right: Corre¬ 
sponding regions (boxes) in the HAADF-STEM image. 


Fluctuation electron microscopy (FEM) 


FEM data were obtained from energy-filtered 
nanobeam diffraction. The intensity fluctuations rep¬ 
resented in the diffraction data were analyzed using 
the annular mean of variance image ^vimage(^) [27]. 
For this purpose, the individual nanobeam diffraction 
patterns (NBDP) of the region of interest were averaged 


HAADF collection semi-angle 
■ 145.8 mrad (48 mm) 

• 87.5 mrad (77 mm) 

51.3 mrad (130 mm) 
t 26.4 mrad (245 mm) 
16.9 mrad (380 mm) 

< 13.4 mrad (480 mm) 



400 

Position / [nm] 


FIG. 5. Dependence of the HAADF intensity on the 

collection semiangle of the HAADF detector. Dark-corrected 
and normalized HAADF intensity line profiles measured for 
different collection semiangles of the HAADF detector along 
a wedge-shaped TEM specimen. The camera lengths are put 
in parentheses. 


pixel-by-pixel to a mean and a variance image, which 
then were converted into annular projections using 
the PASAD (Profile Analysis of Selected Area Diffrac¬ 
tion) tools [39] plugin for Digital Micrograph (Gatan). 
PASAD tools ensure a precise identification of the center 
of the patterns to perform the annular integration. The 
final variance ^vimage(^) was calculated by dividing the 
variance annular average by the mean square. The error 
was estimated using the standard deviation calculated 
from several subsets of patterns. 











































